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Abstract 

The homogenization of a composite material comprising three isotropic dielectric materials was investigated. 
The component materials were randomly distributed as spherical particles, with the particles of two of 
the component materials being coupled to form dimers. The Bruggeman and Maxwell Garnett formalisms 
were developed to estimate the permittivity dyadic of the homogenized composite material (HCM), un¬ 
der the quasi-electrostatic approximation. Both randomly oriented and identically oriented dimers were 
accommodated; in the former case the HCM is isotropic, whereas in the latter case the HCM is uniaxial. 
Representative numerical results for composite materials containing dielectric-dielectric dimers demonstrate 
close agreement between the estimates delivered by the Bruggeman and Maxwell Garnett formalisms. For 
composite materials containing metal-dielectric dimers with moderate degrees of dissipation, the estimates 
of the two formalisms are in broad agreement, provided that the dimer volume fractions are relatively low. 
In general, the effects of intradimer coupling on the estimates of the HCM permittivity are relatively modest 
but not insignificant, these effects being exacerbated by anisotropy when all dimers are identically oriented. 

Keywords: depolarization; homogenization; Bruggeman formalism; Maxwell Garnett formalism; dimer; 
quasi-electrostatic approximation; polarizability density 


1 Introduction 

Composite materials containing random dispersals of particles or inclusions can exhibit remarkable char¬ 
acteristics if two (or more) of their component materials are coupled together to form dimers (or trimers, 
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etc.). Indeed, nanoengineered composite materials containing dimeric particles are playing an increasingly 
prominent role in the development of new technologies [1]. For example, plasmonic interactions in metallic 
dimers can result in enhanced Raman scattering or molecular fluorescence [2, 3, 4], which may lead to highly 
sensitive optical sensors. Dielectric dimers are also of considerable interest, in the context of Mott insulators 
[5, 6], liquid crystals [7], and interstellar molecular hydrogen formation [8], for examples. 

This study concerns the estimation of the effective constitutive parameters of particulate composite ma¬ 
terials, wherein two of the component materials are jointly present as dimers. Previously this topic has been 
investigated using Mie scattering theory [9], a Brownian-motion formulation [10], numerical methods such 
as the finite element method [2], and quantum-mechanical methods based on density function theory [4, 8]. 
In contrast, the theoretical approaches taken herein are simpler, being based on the well-established homog¬ 
enization formalisms named after Bruggeman and Maxwell Garnett [11, 12]. The Bruggeman formalism is 
a self-consistent one in which all component materials are treated in an identical manner. A rigorous basis 
for the Bruggeman formalism arises from the strong-permittivity-fluctuation theory [13]. The nature of the 
Maxwell Garnett formalism is essentially perturbative—which is only appropriate for dilute dispersals of 
particles in a host material [14, 15]. The theoretical basis for the Maxwell Garnett formalism is bolstered 
by its close association with the Hashin-Shtrikman bounds [16]. 

The plan of this paper is as follows. Relevant details of the homogenization formalisms are provided 
in §2. This is followed by the development of appropriate expressions for the electric dipole moments and 
polarizability densities in §3. The homogenization formalisms themselves are set up for composite materials 
containing identically oriented dimers in §4 and randomly oriented dimers in §5. Representative numerical 
results are presented in §6 for both dielectric-dielectric dimers and metal-dielectric dimers. Lastly, a brief 
discussion of the theory and numerical results is provided in §7. As regards notation: the permittivity of 
free space is written as Eq; Cq denotes the speed of light in free space; w is the angular frequency; and 
i = \/^. Vectors are underlined (with the' symbol denoting unit vectors), whereas dyadics [17, 18] are 
double underlined. 


2 Homogenization preliminaries 

The homogenization of three isotropic component materials, labeled as ‘a’, ‘b’, and ‘c’, of a particulate 
composite material is investigated in the following sections. The three component materials are characterized 
by the permittivity scalars Eq, Ef,, and £c- 

Component materials ‘a’ and ‘b’ are dispersed as dimers of electrically small spheres. For simplicity, 
the average radiuses of the spheres belonging to these two component materials are taken to be the same, 
namely s. Each sphere of component material ‘a’ is chemically linked to exactly one sphere of component 
material ‘b’, thus forming a dimer. The distance between the centers of the two spheres making up a dimer 
is d > 2s. In each dimer, the location of the center of the ‘b’ sphere relative to the center of the ‘a’ sphere 
is given by dd, where the unit vector 

d = {ux cos (j) + Uy sin(j)) sin 6 + cos 9, dG[0,7r], ^G[0,27r). (1) 

A heterodimer is specihed by Eq E(,, whereas a homodimer is specified by £a = £b- The dimers are randomly 
distributed. 

When applying the Bruggeman formalism, component material ‘c’ is regarded as a random dispersal of 
electrically small spheres, which we take to have an average radius s. No particular topology is assigned to 
component material ‘c’. in the Maxwell Garnett formalism. The volume fractions of component materials 
‘a’, ‘b’, and ‘c’ are /a, fb, and fc, respectively, with fa = fb and fa + fb + fc = 1- 

Provided that the largest relevant wavelength is much longer than the linear dimensions of the spheres 
and the dimers, the mixture of component materials ‘a’, ‘b’, and ‘c’ may be regarded as being effectively 
homogeneous. The constitutive parameters of the resulting homogenized composite material (HCM) may be 
estimated using a homogenization formalism. Two distinct cases are investigated theoretically: identically 
oriented dimers are the subject of §4, while randomly oriented dimers are treated in §5. 
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All dimers are assumed to be oriented in the same direction in §4; i.e., d is fixed for all dimers. Conse¬ 
quently, the HCM is a uniaxial dielectric material characterized by a permittivity dyadic of the form 

= HCM ^ ^HCM (I - ii) + > (2) 

where / is the identity dyadic. The estimate of (with components and eI/cm) yielded by the 

Bruggeman formalism is written as (with components and e^r) estimate yielded by the 

Maxwell Garnett formalism is written as (with components s^q and £^< 3 ). 

As the dimers are taken to have no preferred orientation in §5, the corresponding HCM is an isotropic 
dielectric material characterized by the scalar permittivity £hcm- The estimate of Shcm yielded by the 
Bruggeman formalism is written as esn and the estimate yielded by the Maxwell Garnett formalism is 
written as £mg- 


3 Electric dipole moments and polarizability densities 

A quasi-electrostatic viewpoint is adopted in both homogenization formalisms, whereby each pair of electri¬ 
cally small spheres comprising a dimer is modeled as a pair of electric dipole moments p^ and p separated 
by the distance d in the direction of d. In the case of the Bruggeman formalism, each electrically small 
sphere of component ‘c’ is modeled as a electric dipole moment p No magnetic dipole moments analogous 
P a b c present in the quasi-electrostatic regime. 

We need to consider two separate cases. In §3.1 the electrically small spheres are immersed in a homoge¬ 
neous uniaxial dielectric material, while in §3.2 the electrically small spheres are immersed in a homogeneous 
isotropic dielectric material. 


3.1 Uniaxial dielectric immersion material 


Let us consider a single dimer immersed in a homogeneous dielectric material characterized by the permit¬ 
tivity dyadic In this section, £^^^ characterizes a uniaxial material whose distinguished axis is parallel 

to d ; i.e., 

= Eout (l-ii) +eLtii- (3) 

The dimer is taken to be centered at the origin r = 0. 

A quasi-electrostatic field, originating from a distant source, is incident on this dimer. The incident 
electric field phasor is denoted by MindL)- In response to E_inciX-)^ electric dipole moments p^ and p^ are 
induced at r ^ = — (d/2) d and r ^ = (d/2) d, respectively. These electric dipole moments are given by [19] 




(4) 


where 



iujD 

out 



(5) 


is the polarizability density dyadic of an isolated sphere of material ‘£’, immersed in the material labeled 

‘out’, with being the corresponding depolarization dyadic [18]. 

The 3x3 dyadics D and a possess the same symmetric form as £ , i.e., 

■’ =out =e/out ■' =Ollt’ ’ 


Rout = (l - ii) + T>Lt M 

^e/out = ^^/out {L-dd)+al/outii ^ 
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The components of are given as [20] 


with 


^■^( 7 ) = ^ 
27 

7 - 1 


DLt = 

7 

L^in) 


£)ll — 

^out 

1 

lII(7) 


1 

tanh” 

U./i 


( 7 ) 


1 - 


: tanh’ 


-'(yrv^) 


VI - 7 

being dimensionless scalar functions of the dimensionless parameter 7 = eoui/^olit- The components of a 

^out ^out) 


( 8 ) 


— tj out 


are given as 


out 


a 


tj out 


eout+7T-^(7) 

p-L fr.-J ^ 

^out I ^outJ 

+ TII(7) (e^-eL) 


^out 


(9) 


The electric field phasor exciting the sphere of material ‘a’ is not merely KinciLa)! instead, 


EexciHa) = ( 10 ) 

where E_^sca represents the electric field scattered by the sphere of material ‘b’. Likewise, there are two 
contributions to ^exciLb)'^ i-e., 

KaxciZb) = Kincij-b) + E^sca(l-b) 7 (H) 

where represents the electric field scattered by the sphere of material ‘a’. In the quasi-electrostatic 
regime, E_inc{T-a) — KindT-b) — ^mc(Q) and the scattered field phasors are given by [17, Sec. 10.5] 


E^scaiZb) 

E\,a{ra) 




47rei-„td3 


2dd — 7 (^L — d^ 

" 

• P 

— a 


> 

2dd — 7 — d^ 

’Pb 

/ 


The combination of Eqs. (4), (10), (11), and (12) yields 



r ^ -v / -v 

4 o \ 

Ea ^^=alout ' 

2dd — 7 fX~ ddj 

'Eb=i^^ ^a,out'^^nM I 


^ ^ / - A 1 

4 I 

Eb-’^^^b/out • 

‘^dd--fU-dd] 

• Ea = 3 ^® ^b/out • ^™c(0) J 


wherein the parameter 

- 

The pair of linear Eqs. (13) deliver the electric dipole moment 

•^mc(Q)) ^e{a,&}, 


( 12 ) 


(13) 


(14) 


(15) 
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with the 3x3 dyadic function 


= 1 ! out 


( 1 - (l- ddj + (l - ddj 

■ [{o^i/out - ^^lOi^/outO^i/out) {L-M) + {<^l/out + 2^-La!/oui«f/o^t) H 


(16) 


The dyadic ^ —e/out polarizability density dyadic of a monomer (sphere) of material as 

a constituent of an isolated dimer immersed in the material labeled ‘out’. The sum 


= dimerjout =ajout =b/out 


(17) 


may be regarded as the polarizability density dyadic of the dimer, and the electric dipole moment 

Edimer ~ Ea~^ Eb ~ 3’’"'® =dimer/out ’ —inci^) 


(18) 


characterizes the quasi-static scattering response of the dimer. Although a uniaxial object [21], the dimer 
is different from a rod or a needle because the volume entering the right side of Eq. (18) is that of a sphere 
but not of a cylinder. 

Next, let us turn to the electrically small sphere of material ‘c’ immersed in a uniaxial dielectric material 
characterized by the permittivity dyadic given in Eq. (3). The sphere is centered at the origin r = 0. 
Suppose that the sphere is illuminated by a quasi-electrostatic field E_inciX-)- The induced electric dipole 
moment is given by [19] 


Ec = • ^-c(0), 


(19) 


where the polarizability density dyadic a ^ ^ is defined per Eq. (5) but with ei therein replaced by Sc-, and 
the components of are written as and per Eq. ( 6 ) 2 - 


3.2 Isotropic dielectric immersion material 


Suppose that the immersion material is isotropic, i.e., = £ outL- Then, the derivations in §3.1 simply 

considerably. In particular, the depolarization dyadic reduces to (SiuiSout) ^ 7, while the polarizability 
density dyadic reduces to a^/outLi where the polarizability density scalar 


^£out ^ 


^out 

££ + 2eoiit J 


£ G {a, b} . 


Consequently, we get 


^ / * 
= ajout 


= bjout 


jout^bjout'} ^a/out^b/out) XX 

O-ab/otii) (l-ii) + (1 + id ^aj out 

ajout^bjou^ XX^ a j out£^h j out) di 


( 1 - cr^a, 
1 - 


-1 ^ 


1 — a‘‘a 


(1 ^£^ajout) dd^ T (t ~t“ ^£££^ajout) 


d d 


^bj out 


( 20 ) 


( 21 ) 


wherein the parameter 


36'oui 


( 22 ) 


These simple expressions are useful when the homogenization of a composite material containing randomly 
oriented dimers is considered. 
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The special case of homodimers is noteworthy. Here £a = £b and thus aa/out = ctbjout- Hence, = 


a, , with 

= bjout 


ia/oui(^) = ^ (= " -"-0 + 

and the dimensionless scalar parameter 


(1 — d) + (1 + 25') dd 


^aj out 1 


G = 


^out \ '5 


(23) 


(24) 


^ “1“ ‘^^out J 

The electric dipole moment of an electrically small sphere of material ‘c’ is given as in Eq. (19) with 
^c/out = ^c/outL where 

Q-c/out = ^Sout ( I”*"* ) ■ (25) 

\ I ^^out J 

4 Identically oriented dimers 

If all dimers have the same orientation then the HCM is a uniaxial dielectric material [17, 18] with its 
distinguished axis parallel to the fixed unit vector d. That is, has the form given in Eq. (2). 

4.1 Bruggeman formalism 

Particles of all component materials are assumed as being immersed in the HCM itself, in the Bruggeman 
formalism [11, 12]. Thus, the expressions presented in §3.1 are appropriate here with the subscript ‘Hr’ 
replacing the subscript ‘out’. The Bruggeman formalism rests on upon the assumption that the electric 
dipole moments arising from the electrically small spheres of the component materials, weighted by volume 
fraction, sum to the zero vector [22]; i.e., 


f T) 

^ ^ dimer 


+ fcP, = fa{E^+P,)+fcP^ = 0. 
Upon combining Eqs. (15)-(19) with Eq. (26), the dyadic equation 
= “ f‘^=dimer/Br~^ ^‘^=c/Br 

~ f°'^a/Br~^ ■^’^^b/Br^ ^^=c/Br 


(26) 

(27) 


/a ( 1 - <^^l^a^/Brat/Br) (l “ M) + (l “ H 

* (*^a/Br ^bjBr Br^bjBr^ (= dd^ 

+ (' 


Q!"/d +o}},in +4crQ;^^D aLo I dd 

a/Br b/Br a/Br h/Br J - 


+ /c 


^c/Br 


(l - ddj 


+ al/B^dd 


(28) 


emerges. 

Due to the uniaxial symmetry, Eq. (28) represents two coupled scalar equations with and as the 
two unknowns, which must be obtained by numerical methods. The following Jacobi scheme may be used for 
this purpose [24]. First, let us notice that t/Br’ defined in Eq. (16) (with the subscript ‘Hr’ replacing 
the subscript ‘out’), may be written as 


/n = ^ /n 

= ajBr ^=0 =ajBr 

a./n = M_ -a 

= bfBr ^=a =bfBr 


(29) 
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wherein the 3x3 dyadic 




Add 


' ^ /n 
= al Br 


^b/Br} 


|i+(T (^L-ddj +2dd i£{a,b}. 


(30) 


Second, notice that as defined in Eq. (5) (with the subscript ‘Br’ replacing the subscript ‘out’), 

may be written as 

wherein the 3x3 dyadic 


= i/Br 


{etl - g 


■ P 

Br =t^ 


i G {a, b, c} , 


(31) 


-1 


(32) 


g,= [L+^U:^Br' (^4-iB. 

Hence, after using Eqs. (29) and (31), Eq. (27) may expressed as 

/a[M, • (eaL-gs^) -^a+Ma* {^bL-gBr) * £ J + A (^c " £ ^ J 'E,=g- (33) 

The Bruggeman estimate may be extracted from Eq. (33) by the iterative scheme represented by 

= Br ~ ^ {=Br} ’ 

where the action of the dyadic operator T is given by 


(34) 


r 


{sb,} = [/-(Mr £.+M. •£») + /=£< 


£(£.Mr£. + s*Mr &)+£="£, 


(35) 


4.2 Maxwell Garnett formalism 


Particles of component materials ‘a’ and ‘b’ are viewed as immersed in component material ‘c’, in the Maxwell 
Garnett formalism [11]. Thus, the expressions presented in §3.2 are appropriate here with the subscript ‘c’ 
replacing the subscript ‘out’. The electric dipole moments of spheres of component material ‘c’ are not 
relevant to this formalism, the HCM essentially arising as a perturbation of the component material ‘c’ by 
the addition of a relatively small amount of component materials ‘a’ and ‘b’. Consequently, results of the 
Maxwell Garnett formalism, as applied here, are strictly appropriate only for fa ~ 0.15. 

The Maxwell Garnett estimate of is obtained explicitly as [23] 


= MG 


^cL~\- fa 



A — — 01 , 

= 3er =“/= 


- a 


h/c 



2/a - 

3£a =^/- 


(36) 


with the polarizability density dyadics as given in Eqs. (21). 


5 Randomly oriented dimers 

If the dimers are randomly oriented, the HCM is isotropic. Accordingly, the electrically small spheres of 
component materials ‘a’ and ‘b’ (and ‘c’ in the case of the Bruggeman formalism) should be regarded in 
this case as being immersed in an isotropic dielectric material and the expressions presented in §3.2 are 
appropriate. 
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Orientationally averaged electric dipole moments are defined as 


^ n2iTZ pTT 

{p ) = — / p sine dedcj), iG{a,b,c}. 


( 37 ) 


277 pTT 

= 0 Jl9^0“ 

The orientationally averaged electric dipole moments for component materials ‘a’ and ‘b’ may be expressed 

47r 


as 


(P,) = ys" sin0d0d</>) •^„,(0), iG{a,b}. 


(38) 


Herein the quantity in parenthesis represents the orientational average of the polarizability density dyadic 


a„, ; for later use, this is written as 

=ifout 




ave T 
Ijout ^ ’ 


i £ {a, b} . 


Let us also note that 


(Ed^mer'^ = (Pn) + iP,) 


tlb' 


'' = dimer/out^ ^^a/out^ ^^b/out^ 
ave _ ~ave , -ave 

dimer j out afout ' bjout 


> . 


Upon evaluating the integrals on the right side of Eq. (38), the following result is delivered: 


9 ^*^ ^ajout 


(P^) = -^^S^Oib/out 


2 (l ^Ca^b/out') 


1 H" ‘^^^bjout 


1 (J (^ajout^bjout f d(J OLq^IQ iiiOL}yjQy^i 

2 (l ^ 1 “1“ 2(701 aj out 

1 (7 j out^b I out f ^ajout^bjout 


E^uM) 


E,nM 


. 


(39) 


(40) 


(41) 


In the special case of homodimers, P ^ = P the corresponding orientationally averaged electric dipole 

moment is given as 


^Pj = ^Pb> = o^S^«a/o«t 


d — 1 


(ct + 1) (2cr — 1) 


E.uM). 


(42) 


Since spheres of material ‘c’ have no directional dependency, the orientational average of the associated 
electric dipole moment is simply p itself; i.e.. 


(p) = ■^'^S^OL^/outE^nM)- 


(43) 


5.1 Bruggeman formalism 

In the Bruggeman formalism, particles of all component materials are assumed as being immersed in the 
HCM itself. Thus, the expressions presented in §3.2 can be used here with the subscript ‘out’ replaced by 
the subscript ‘Br’. The Bruggeman formalism dictates that 

= S' (44) 

Upon combining Eqs. (41) and (43) with Eq. (44), the corresponding scalar Bruggeman equation emerges as 


j ^a{Br 

2(l — aab/Br) , l + 2cra{,/Br 

1 3 

O^b/Br 


1 (7 O^djBr^bjBr ^ ^ajBr^bjB 

2 (l Br^ 1 f “^^^a j Br 

r 

3 

1 (7 CV-dj^j'CXbfBr f d(7 ^ajBr^bjBr 
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( 45 ) 

























After using Eqs. (14), (20), and (25) to substitute for cr, Ua^b/Br, and ac/Br, respectively, Eq. (45) may 
be recast as a quintic polynomial in Ebt- (In the case of homodimers, this polynomial reduces to a cubic 
polynomial in Ebt)- A Jacobi numerical scheme [24] can be used to extract Ebt from Eq. (45). That is, the 
solution may be found using the iterative scheme represented by 


EBt — ^ : 

where the action of the scalar operator S is given by 

TTlbSa , '^a 


+ 


fc £c 


S{SBr} = 


£a + ‘^SBr £b + “^.EBr Ec + ‘2:£B 


rUb 


+ 


fc 


with the scalar parameters 


£o + ‘^SBr £b + Sesr £c + ‘^^Br 

2(1 — Craa/Br) ^ + ‘^<^0'a/Br 


mb = - 


1 fJ Oi(ijBr^bjBr 1 4fJ bXajBr^bjBr 

2 (l EEOibjBT^ 1 4“ ‘^EbO^bjBr 


1 E^^OtajBr^bjBr 1 Br^bj Br 


. 


(46) 


(47) 


(48) 


5.2 Maxwell Garnett formalism 


In the Maxwell Garnett formalism, particles of component materials ‘a’ and ‘b’ are viewed as immersed in 
component material ‘c’. Accordingly, here the expressions presented in §3.2 are used with the subscript ‘c’ 
replacing the subscript ‘out’. The Maxwell Garnett estimate of Ehcm is given explicitly by 


£mg — £c s 1 + 3/a 


■b/c 


3ea-2/ad“7, 3£c-2fa 


a 


b/c 


where 


a. 



- ave _ 

^b/c ~ '^^b/c 


2 (l - aab/c) 

1 fJ O^ajc^bjc 

2 (1 - aag/c) 

1 fJ O^djcCXbjc 


1 + 2aab/c 

I 40" E^ajc^bjc 

1 + 2(jaa/c 

I 4fJ OiajcG-bjc 


. 


(49) 


(50) 


6 Numerical results 

Let us now present representative numerical results based on the theoretical results established in §2-§5. 
Although the range fa G [0, 0.5] may appear appropriate at first glance, the maximum value of fa must be 
less than 0.5. This is because no sphere of material ‘c’ must be allowed to occupy the region between the 
two spheres constituting a dimer. Nevertheless, we have provided the Bruggeman estimates for fa G [0,0.5], 
because the upper limit of fa will have to be decided experimentally for a specific composite material. The 
Maxwell Garnett formalism is appropriate only for dilute composite materials, and we have restricted the 
presentation of the Maxwell Garnett estimates to fa G [0,0.15]. 

6.1 Dielectric—dielectric dimers 

Suppose, first, that both component materials ‘a’ and ‘b’ are nondissipative dielectric materials, specified 
by the permittivities £a = 2£a and £b G (£ro,10£rQ). The dielectric-dielectric dimers which arise from the 
combination of component materials ‘a’ and ‘b’ are randomly dispersed along with component material ‘c’ 
specified by the permittivity Ec = (14 + ii) Eg. 
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6.1.1 Randomly oriented dimers 

The real and imaginary parts of SmgI^o and esrl^o are plotted against Eb/^o and fa in Fig. 1, for the case 
where the dimers are randomly oriented and d = 2s. The real part of Emg decreases approximately linearly 
as fa increases, with its rate of decrease being greatest at the lowest values of Eb- Furthermore, the real part 
of Emg increases approximately linearly as Eb increases, with its rate of increase being greatest at the largest 
values of fa- The imaginary part of Emg decreases approximately linearly as fa increases, and this trend is 
almost independent of the value of Eb- For the range 0 < /a ~ 0.15, the real and imaginary parts of Esr are 
very similar, both qualitatively and quantitatively, to the real and imaginary parts of emg- For fa ~ 0.15, 
both the real and imaginary parts of Esr exhibit a more nonlinear dependency on fa than they do at lower 
values of fa- 

The issue of the influence of the dimer separation distance upon Ehcm is addressed via Fig. 2, wherein 
the real and imaginary parts of EmgI^o and EbtIeq are plotted against d/s for Eb = Eq (green, solid curves), 
2eo (blue, dashed curves), and Se,, (red, broken dashed curves), when fa = 0.15. The influence of d on the 
real and imaginary parts of Ehgm, for both the Bruggeman and Maxwell Garnett estimates, decays rapidly 
as d increases. Indeed, both estimates of Ehcm are practically independent of d for d > 4s, these estimates 
being essentially the same as those that would be obtained through the homogenization of three independent 
component materials ‘a’, ‘b’, and ‘c’ with no dimeric interaction between the spheres of component materials 
‘a’ and ‘b’. The change in the real part of Emg as d increases from zero to 4s is approximately 0.2%, whereas 
the corresponding change in the imaginary part of Emg is approximately 0.4%. The Bruggeman estimates 
of Ehgm are somewhat more sensitive to changes in d, the real part of Ebt changing by approximately 1.6% 
whereas the imaginary part of Ebt changing by approximately 2.8% as d increases from zero to 4s. 

6.1.2 Identically oriented dimers 

Qualitatively, the Bruggeman and Maxwell Garnett estimates of e^^qj^ and Ej^Qj^ for composite materials 
containing identically oriented dimers have dependencies similar to the estimates of Ehcm presented in 
Figs. 1 and 2 for composite materials containing randomly oriented dimers. This becomes evident from the 
plots of the averages (^e^mq +£mg) with respect to Eb/Eg and fa in Fig. 3. To 

the naked eye, the plots in Figs. 1 and 3 are almost indistinguishable. 

However, there are significant quantitative differences between the estimates of e\^qj^ and E'^qm- 
Fig. 4, the real and imaginary parts of the differences (^e^mq — Ej^Q^ /eq and /eq are plotted 

against Eb/Eg and fa, when d = 2s. Both the real and imaginary parts of the difference E^j^^Q — E-j^Q increase 
approximately linearly as fa increases and decrease approximately linearly as Eb increases. Thus, the greatest 
degree of anisotropy is predicted by the Maxwell Garnett formalism when fa is largest and Eb is smallest. 

In the range 0 < /a ~ 0.15, both the real and imaginary parts of the difference — e^^, are qualitatively 
similar to the corresponding real and imaginary parts of e^Ijq — However, at larger values of fa, both 

the real and imaginary parts of e^,^ —e;^^ exhibit strong nonlinear dependencies on fa. The greatest degree of 
anisotropy is predicted by the Bruggeman formalism to be in the vicinity of fa ~ 0.2 with Eb = Eq. Broadly, 
over the parameter ranges considered, the degree of anisotropy estimated by the Maxwell Garnett formalism 
is slightly larger than that estimated by the Bruggeman formalism. 

The degree of anisotropy exhibited by the HCM, as estimated by the Bruggeman and Maxwell Garnett 
formalisms, decays rapidly as the separation distance d in the dimer increases. This is illustrated in Fig. 5 
wherein the real and imaginary parts of the differences ^£^<5 — Emg) /^o and ^£ 5 ^ — E^r) /^o are plotted 
against d/s for Eb = £0 (green, solid curves), 2 £o (blue, dashed curves), and Seq (red, broken dashed curves). 
Here fa = 0.15. The graphs for e^^ — Eg^ and eI^q — e/^q in Fig. 5 are qualitatively similar, with the 
Maxwell Garnett estimates being slightly larger than the Bruggeman estimates at each value of d and Eb. 
The degree of anisotropy, as estimated by both formalisms, falls most rapidly for the smallest value of £&. 
Furthermore, the degree of anisotropy, as estimated by both formalisms, vanishes almost entirely at d = 4s. 
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6.2 Metal—dielectric dimers 


Next, suppose that component material ‘a’ is a metal. For definiteness, this metal is taken to be silver as 
characterized by the size-dependent permittivity [25, 26] 


£Ag{s) — So 


1 - 


lAg 


3vf 

4s 


(51) 


Herein, vp = 1-4 x 10®m is the electron speed at the Fermi surface, jAg = 10^"^ s“^ is the relaxation 
rate, and tOp = 1.38 x 10^® rad is the plasma frequency. The angular frequency is a; = 27rco/Ao, with 
the free-space wavelength chosen to be Ao = 650 nm. As in §6.1, component material ‘b’ is a nondissipative 
dielectric material specified by the permittivity e& G (so, IOeq)- The metal-dielectric dimers which arise from 
the combination of component materials ‘a’ and ‘b’ are randomly mixed with component material ‘c’ which 
is specified by the permittivity Sc = (14 -|- 4i) Sg for all results presented here. 

6.2.1 Randomly oriented dimers 

For composite materials containing randomly oriented dimers, the real and imaginary parts of Emg/eo a-nd 
SBrIso are plotted in Fig. 6 against Sb/So and fa for d = 2s and s = 5 nm. Thus, Sa = (—21.4 -|- 2.4i)eo 
by virtue of Eq. (51). The graphs of the real parts of Smg and ssr in Fig. 6 are qualitatively similar to 
the corresponding graphs in Fig. 1 for dielectric-dielectric dimers. In contrast, graphs of the imaginary 
parts of Smg and Sbt in Fig. 6 are rather different to the corresponding graphs in Fig. 1, both qualitatively 
and quantitatively. On the whole, the imaginary parts of Smg and Sbt are substantially larger in Fig. 6 
than they are in Fig. 1. Furthermore, the imaginary parts of Smg and Sbt in Fig. 6 are substantially more 
nonlinear with respect to increasing fa than are the corresponding quantities in Fig. 1. Both the real and 
imaginary parts of Smg and SBr are qualitatively similar in Fig. 6 in the range 0 < /a ~ 0.15. However, 
across this range, the quantitative differences between the estimates Smg and Sbv are substantially larger 
than the corresponding differences presented in Fig. 1, and these differences between smg and SBr grow as 
fa increases from zero. 

Qualitatively, the influence of the dimer separation distance d upon Shgm for the metal-dielectric dimer 
HCM is similar to that for the dielectric-dielectric dimer HCM. This may be appreciated by comparing 
Fig. 2 with Fig. 7. In Fig. 7, the real and imaginary parts of Smg/so and SbvIso are plotted against d/s 
for Sb = Sq (green, solid curves), 2eo (blue, dashed curves), and SsTo (red, broken dashed curves), for the 
metal-dielectric dimer HCM, with fa = 0.02 and s = 5 nm. As is the case for dielectric-dielectric dimers 
in Fig. 2, the estimates of Shgm for metal-dielectric dimers in Fig. 7 are practically independent of d for 
d > 4s. The magnitudes of the relative changes in Shgm as d increases from zero to 4s in Fig. 7 are similar 
to those in Fig. 2, with the Bruggeman estimates being somewhat more sensitive than the Maxwell Garnett 
estimates to changes in d. 

The effects of the size of the metal particles which make up component material ‘a’ are delineated in 
Fig. 8. Herein the real and imaginary parts of Emg/sq and SBr/sg are plotted against s for Sb = £0 (green, 
solid curves), 2£o (blue, dashed curves), and 3£o (red, broken dashed curves), with fa = 0.02 and d = 2s. 
The influence of s on the real and imaginary parts of Shgm, for both the Bruggeman and Maxwell Garnett 
estimates, decays rapidly as s increases. Indeed, the estimates of Shcm vary little as s increases beyond 
20 nm. The change in the real part of Smg as s increases from 5 to 20 nm is approximately 2.6%, whereas 
the corresponding change in the imaginary part of Smg is approximately 1.0%. The real part of Sbt changes 
by approximately 1.2% whereas the imaginary part of Sbt changes by approximately 1.8%, as s increases 
from 5 to 20 nm. Most strikingly, the real and imaginary parts of Smg, as well as the real part of Sbt, 
uniformly decrease as s increases from 5 to 20 nm whilst the imaginary part of eBr uniformly increases as s 
increases. 
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6.2.2 Identically oriented dimers 

The estimates of and yielded by the Bruggeman and Maxwell Garnett formalisms for composite 

materials containing identically oriented dimers exhibit characteristics that are qualitatively similar to those 
displayed in Figs. 6-8 by the corresponding estimates of £hcm for randomly oriented dimers. However, 
significant quantitative differences arise between the estimates of ^hcm ■ imaginary 

parts of the differences — ^mg) /^o and /£„ are plotted in Fig. 9 against Eb/so and fa, 

for d = 2s and s = 5 nm. Across the range 0 < /„ £ 0.15, the differences ef/GM ~ ^hcm estimated by 
the Maxwell Garnett and Bruggeman formalisms are qualitatively similar. However, there are quantitative 
differences between £ 3 ,^ ^ad £\^q ~£mg^ and these increase in magnitude as fa increases. The greatest 
degree of anisotropy is estimated by the Bruggeman formalism to exist when both fa and £b are maximum. 
In contrast, the greatest degree of anisotropy is estimated by the Maxwell Garnett formalism to exist when 
fa is maximum but Eb is minimum. 

The influences of the dimer separation distance d and the metal sphere radius s upon the anisotropy of 
the HGM are delineated in Fig. 10. Here the real and imaginary parts of the differences ^£|^g ~ ^mg) /^o 

and — Eg^'^ /eq are plotted against s (in nm) and d/s, for fa = 0.02 and d = 2s. While the real part 

of e\jq — e^q decreases uniformly as d increases from 2s to 4s, this quantity varies very little as s increases 
from 5 nm to 20 nm. The imaginary part of £|^g ~ ^mg increases uniformly as d increases from 2s to 4s; in 
contrast, the imaginary part of e\jq — e/^q varies only marginally as s increases from 5 nm to 20 nm. The 
graphs for the real and imaginary parts of £^^ — Eg^ are both qualitatively and quantitatively similar to the 
corresponding graphs for e\^q — e-^q . 

7 Discussion 

The Bruggeman and Maxwell Garnett formalisms have been established in the preceding sections for the 
homogenization of composite materials containing randomly oriented and identically oriented dimers. The 
representative numerical results presented in §6.1 for the case of dielectric-dielectric dimers demonstrate 
close agreement between the estimates of the HGM constitutive parameters delivered by the Bruggeman and 
Maxwell Garnett formalisms for both randomly oriented and identically oriented dimers. The Bruggeman 
formalism is advantageous over the Maxwell Garnett formalism insofar as the former is appropriate for 
arbitrary dimer volume fractions whereas the latter is appropriate only for low dimer volume fractions. On 
the other hand, the Maxwell Garnett formalism is relatively straightforward to implement numerically as 
its estimates are provided as explicit formulas, in contrast to the Bruggman formalism whose numerical 
implementation typically involves the iterative extraction of estimates from implicit formulas. 

The case of metal-dielectric dimers should generally be approached with caution. If attention is restricted 
to parameter regimes involving low dimer volume fractions and moderate degrees of dissipation then, as 
demonstrated in §6.2, the Bruggeman and Maxwell Garnett formalisms deliver estimates of the HGM con¬ 
stitutive parameters which are in broad agreement. However, at higher dimer volume fractions, substantial 
qualitative and quantitative differences emerge between the estimates provided by the two formalisms, and 
these differences are exacerbated by anisotropy in the case of identically oriented dimers. 

The difficulties that arise for metal-dielectric dimers at larger values of the dimer volume fraction es¬ 
sentially stem from the fact the real part of Ea is negative while the real parts of Eb,c are positive. In the 
absence of substantial degrees of dissipation, homogenization for such parameter regimes can be problem¬ 
atic for conventional formalisms, especially at mid-range values of volume fractions [27, 28, 29, 30]. For 
examples, in these parameter regimes the Bruggeman estimates may violate the Hashin-Shtrikman bounds, 
and the Maxwell Garnett estimates may exhibit very large resonances (with respect to varying volume frac¬ 
tion). These issues effect both passive and active HGMs [31], for both forward and inverse homogenization 
formalisms [32], and also impose limitations on the Bergman-Milton bounds [33]. 
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The mathematical origin of these problematic parameter regimes may be appreciated most readily by 
considering the simplest case, namely that of randomly oriented homodimers. The corresponding expression 
for the orientationally averaged electric dipole moment is provided in Eq. (42). The polarizability scalar 
Ota/out therein becomes infinitely large in magnitude in the limit £a —t —2£out- This eventuality — which is 
sometimes referred to as a Frohlich mode [25] — may result in singular or highly resonant behavior in the 
estimates of HCM permittivity. There is further scope for singular behavior which is solely attributable to the 
dimer interaction: the denominator term in Eq. (42) is null valued at d £ {~1,1/2}. By using the definition of 
a provided in Eq. (24) with d = 2s, these singularities arise in the limits £a —>■ —(17/7)eout and £a —>■ S^out, 
respectively. For strictly nondissipative materials, the possibility of Salsout S {—2,—(17/7),—3} can only 
arise if either £a£b < 0 or £a£c < 0. Thus, by extrapolation, it may anticipated that regimes in which the 
real part of Eq is negative while the real parts of at least one of £& or Sc is positive may well be problematic. 
However, as demonstrated in §6.2, provided that only low dimer volume fractions are considered and there 
is a moderate degree of dissipation, the problems of singular or highly resonant behavior may not arise. 

The numerical results in §6 reveal that the effects of intradimer coupling decay rapidly as d increases. 
Indeed, for d > As these effects are generally negligible and the permittivity dyadic of the HCM is practically 
the same as that which would arise in the case where component materials ‘a’ and ‘b’ were not coupled at 
all. By comparing the numerical results at d = 2s with those at d = 4s in Figs. 2, 5, 7, and 10, it may be 
deduced that intradimer coupling generally has relatively modest but not insignificant effects on the HCM 
parameter estimates delivered by the Bruggeman and Maxwell Garnett formalisms, and that these effects 
are exacerbated by anisotropy in the case of identically oriented dimers. 

Acknowledgement. AL thanks the Charles Godfrey Binder Endowment at Penn State for ongoing 
support of his research activities. 
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Figure 1: The real and imaginary parts of EmgI^q and EbtIeo plotted against Eb/Eg and fa for the case where 
the electrically small spheres of component materials ‘a’ and ‘b’ combine to form dielectric-dielectric dimers 
with Ea = 2^0 and Eb £ (1,10) Eg- Component material ‘c’ is specified by Eg = (14 + 4i) eo. The dimers are 
randomly oriented and d = 2s. 
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Figure 2: As Fig. 1 except that the real and imaginary parts of Emg/eo and esr/eo are plotted against d/s 
for Ea = 2^0 and et = Eo (green, solid curves), 2eo (blue, dashed curves), and Seq (red, broken dashed curves). 
Here fa = 0.15. 
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Figure 3: As Fig. 1 except that the dimers are identically oriented and the real and imaginary parts of the 
averages i£\jQ + £mg) + ^Br) are plotted. 
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Figure 4: As Fig. 1 except that the dimers are identically oriented and the real and imaginary parts of the 
differences — Smc} /^o and /eq are plotted. 
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Figure 5: As Fig. 2 except that the dimers are identically oriented the real and imaginary parts of the 
differences (s^mq — Smg] /^o and ) /eg are plotted. 
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Figure 6: The real and imaginary parts of EmgI^q and EbtIeo plotted against £b/eo and fa for the case where 
the electrically small spheres of component materials ‘a’ and ‘b’ combine to form metal-dielectric dimers 
with Ea = (—21.4+ 2.4i)eo (i.e., Ea = eAg(5nm) for Aq = 650 nm) and Eb £ (1,10) Eq- Component material 
‘c’ is specified by {ec/eq) = 14 + 4z. The dimers are randomly oriented and d = 2s. 
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Figure 7: As Fig. 6 except that the real and imaginary parts of Emg/^o and esrl^o plotted against d/ 
Eb = £o (green, solid curves), 2£o (blue, dashed curves), and Sso (red, broken dashed curves). Here fa = 
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Figure 8: As Fig. 7 except that the real and imaginary parts of emgI^q and esr/^o are plotted against a 
(nm) with d = 2s. 
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Figure 9: As Fig. 6 except that the dimers are identically oriented and the real and imaginary parts of the 
differences — Smc] /^o and /eq are plotted. 
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Figure 10: As Fig. 9 except that — Emg) /^o (^Sr ~ ^Sr) /^o are plotted against s (nm) and d/s. 

Here £& = 2eo and fa = 0.02. 
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